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Abstract 

We investigated the nanosecond-scale time decay of the blue-green light emitted by nominally 
pure SrTiC>3 following the absorption of an intense picosecond laser pulse generating a high density 
of electron-hole pairs. Two independent components are identified in the fluorescence signal that 
show a different dynamics with varying excitation intensity, and which can be respectively modeled 
as a bimolecular and unimolecolar process. An interpretation of the observed recombination kinetics 
in terms of interacting electron and hole polarons is proposed. 
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I. INTRODUCTION 



Transition-metal oxides with a perovskite-type structure form a fascinating class of mate- 
rials, with an extraordinarily varied physics. Among them, SrTiC>3 (STO) is a prototype ma- 
terial, with a simple cubic perovskite structure at room temperature. Its electrical properties 
appear deceptively simple: it is a band-gap insulator, with a very large dielectric constant due 
to the high polarizability of its ionic lattice. Mobile charges can be however added to STO 
by chemical doping or photon excitation, turning it into a polaronic conductor 1, Q and, at 
very low temperatures, even into a superconductor. [3 ISTO interfaces with other polar oxides 
have recently shown surprising transport properties, % Q in connection with charge transfer 



through the interface and with polaron localization effects. [6| Even in bulk S 



nature of the electron and hole polarons is still fairly unclear. 
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O, the precise 



11 



12| As for 



its optical properties, STO has recently gained renewed attention following the discovery of 
a significant blue-light photoemission from Ar + -irradiated or n-doped STO samples. Q] 



This finding, besides its interest for potentia' 



nescence phenomenology of this material. 



15 



app 



16 



ications, adds to the rather puzzling lumi- 
17, Q] The photoluminescence from a 



material is also a direct probe of the dynamics and recombination processes of its photoex- 
cited charge carriers. 

Motivated by these considerations, in this work we investigated the nanosecond-scale 
time-resolved fluorescence of nominally pure STO, following excitation from an intense pi- 
cosecond laser pulse generating a "plasma" of electrons and holes having initial densities well 
above 10 20 cm~ 3 . To our knowledge, no study probing this regime was reported hitherto, 
for STO or any other perovskite-type oxide. To avoid any possible misunderstandings, we 
should however emphasize that, in this paper, we are using the term "dense plasma" in a 
generic sense, i.e., to indicate a system with a very large number of interacting electrons 
and holes that are, at least initially, freely mobile (although it is possible, in our case, that 
at least a fraction of them will localize a short time after being generated), regardless of 
the strength of their electrostatic interaction and of the clear manifestation of many-body 
effects. It must be also noted that photoluminescence spectroscopy in this regime, with a 
photogenerated carrier density much higher than impurity concentration, is likely probing 
mainly intrinsic properties of the material. 

This article is organized as follows. Section II describes the samples and the experimental 
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procedure. The experimental results and their modeling are reported in Section III. These 
results are then discussed and tentatively interpreted within a microscopic framework in 
Section IV. Section V contains some conclusive statements. 



II. EXPERIMENT 

We investigated five transparent stoichiometric STO crystals produced by four different 
companies (all with specified impurity levels below 150 ppm), in the form of 1 mm thick 
plates. All measurements discussed in this paper were taken at room temperature. In 
our photoluminescence experiments, the excitation was induced by ultraviolet (UV) pulses 
having a wavelength of 355 nm, obtained as third-harmonic of the output of a mode-locked 
Nd:YAG laser. The corresponding excitation photon energy is 3.49 eV, well above the 



indirect band-gap of STO (3.26 eV) and just above its direct-transition edge (3.46 eV).|20l.l21| 
The UV pulses had a duration of 25 ps (full-width at half-maximum), a repetition rate of 
10 Hz, and an energy of — 2 mJ. The UV beam was weakly focused on the sample, with 
a spot size of 1.2±0.1 mm (radius at 1/e 2 of the maximum). The resulting fluence in the 
spot-center was therefore — 90 mJ/cm 2 . Since at a wavelength of 355 nm the optical 
penetration depth is of about 1 fim and the reflectance is about 25%, 20j, |2lJ] the estimated 
peak (spot-center) energy density absorbed in the material surface layer is ~ — 600 J/cm 3 . 
This corresponds to a peak density of photo-generated electron-hole (e — h) pairs of up to 
1.2 xlO 21 cm -3 . During all our experiments, no visible photoinduced damage of the sample 
surfaces was induced and no irreversible variation of the signal with time was seen. 

The luminescence emitted from the sample was collected by a lens system imaging the 
illuminated sample spot onto the detector head, after blocking the (much stronger) elas- 
tic scattering by a long-pass filter with a cutoff wavelength of 375 nm. A typical (time- 
integrated) luminescence spectrum is shown in the inset of Fig. [U It peaks at 425 nm, in the 
blue, corresponding to 2.9 eV of photon energy, with a long tail reaching about 650 nm (1.9 
eV). All samples, independent of the manufacturer, had a roughly comparable fluorescence 
yield, with a maximum sample-to-sample variation of 80%. The initial emission intensity 
was found to be even more stable, with a sample-to-sample variation below 30%. This 
stability is a clear indication that this luminescence is intrinsic, although a dependence on 
intrinsic lattice defects (i.e., not associated with impurities) cannot be excluded. Although 
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FIG. 1: (Color online) Example of photoluminescence signal. Main panel: luminescence decay 
(pulse energy of 550 /xJ). Data are shown as (red) filled circles, the black solid line is the the- 
oretical fit, taking into account the instrumental response (the latter causes some minor signal 
features, such as the small wiggling seen at 10 ns delay). The unimolecular (UD) and bimolecular 
(BD) contributions (see text for a definition) are shown as (blue) dashed and dot-dashed lines, 
respectively. Inset: emission spectrum (the narrow peak seen at 355 nm is due to residual UV 
elastic scattering). 

we have not performed absolute yield measurements, the photoluminescence is clearly weak, 
indicating the dominance of non-radiative energy relaxation. Finally, in all samples the 
spectrum profile and the overall yield were found to be independent of excitation intensity, 
with no significant sign of saturation, up to about 1 mJ of pump pulse energies. We see in- 
stead a significant decrease of the yield above 1 mJ, as shown in Fig. [2j This extremely high 
saturation intensity clearly cannot be ascribed to the filling of impurity levels and therefore 
is a further confirmation of the intrinsic nature of this photoluminescence. This saturation 
might be instead related with the onset of the Auger recombination effect (we will discuss 
this possibility later) or with plasma electromagnetic screening effects (while we may prob- 
ably exclude the contribution of absorption "bleaching", which would increase the optical 
penetration length and therefore reduce the absorbed energy density, but would not reduce 
the overall absorbed energy and hence the luminescence emission). The estimated plasma 
frequency of the photogenerated carriers, assuming free electron masses for both electrons 
and holes, is of about 6 x 10 14 Hz at the highest excitation intensities, corresponding to an 
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FIG. 2: (Color online) Photoluminescence integrated emission intensity /tot versus excitation pulse 
energy E. Square dots are data. The solid line is a linear best-fit in the energy range 0-1 mJ. The 
dashed line is a best-fit based on the typical saturation behavior I tot txE/(l + E/E s ), that yields 
E, = 1.3 mJ. 



electromagnetic wavelength of 500 nm. Therefore it is possible that some partial screen- 
ing of the UV excitation light is taking place at these energies, thus decreasing the sample 
absorption efficiency. In the rest of this paper, our analysis will be mainly focused on the 
non- saturated energy range 0-1 mJ. 

In time-resolved measurements, the luminescence was detected with a photodiode (PD) 
having a rise-time of about 150 ps. In most experiments the entire luminescence spectrum 
was integrated (with a weight given by the PD detection efficiency), but in some cases we 
inserted another filter in the detection line in order to select a portion of the spectrum. The 
PD signal was acquired by a 20 Gsample/s digital oscilloscope having an analog bandwidth of 
5 GHz. With careful data analysis, taking into account the instrumental response function 
(see below), this set-up has a time resolution of about 100 ps. An example of a typical 
measured decay of the luminescence is shown in Fig. [TJ The acquired signal S(t) is actually 
the result of a convolution between the real luminescence decay I(t) and the detection- 
system response function r(t), where t denotes time. The response function was measured 
by detecting the elastic-scattered light of the excitation pulse. Once r(t) and S(t) are known, 
the actual I(t) decay may in principle be calculated by a deconvolution procedure. However, 
in order to avoid the typical problems associated with numerical deconvolutions of data, we 
prefer here to compare the unprocessed raw signals S(t) with the theoretical luminescence 
decays convoluted with r(t). 
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FIG. 3: (Color online) Photoluminescence decay dependence on excitation pulse energy (in the 
non-saturated range 0-1 mJ). Both panels (a) and (b) show the same decay data (circles) referring 
to one of the investigated samples, with the corresponding best-fit curves (black solid lines) as 
emerging from a single global fit (see text). Data and fits are normalized to their maximum. Panel 
(a) is in linear scale and different decays are vertically offset for clarity. Panel (b) is in semi- 
log scale, with no offsetting. The inset in panel (b) is a zoomed-in view on the initial decay. The 
oscillations in the best-fit curves seen in panel (b) result from the convolution with the instrumental 
response function. 

III. RESULTS AND ANALYSIS 

When varying the excitation intensity, the measured decay profile changes markedly, as 
shown in Fig. [31 In particular, the initial decay becomes clearly faster for higher excitation 
intensities (see also the inset in Fig. [3b), while the slower "tail" changes only its amplitude, 
but not its decay rate. This dependence on excitation intensity excludes an interpretation 
as a simple "unimolecular" exponential decay (UD) process and implies that particle corre- 
lations are relevant. A "bimolecular" non-exponential decay (BD) naturally arises from a 
model in which holes and electrons (or the corresponding polarons) have a "partner-capture" 
probability per unit time that is proportional to the density of partners, and this capture 
event is also the rate-limiting step for e — h recombination (i.e., it may correspond either to a 
direct recombination or to a recombination via an intermediate short-lived state). Labeling 
with N(t) the (equal) density of photo-generated electrons and holes, the simplest "pure 
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BD" model is described by the following rate equation 

dN 9 . . 

~dt = «' « 

where 7 is a constant. For excitation pulse energies below the saturation level, the initial 
condition to be associated to this equation is N(0) = all, where U is the excitation energy 
density and a is a constant. We tried this model and found that it provides an excellent 
fit to the energy dependence of the initial faster decay but cannot account for the slower 
tail. An obvious extension of the model is then to introduce also a UD contribution in 
the decay. A unimolecular recombination arises when there is decay channel in which the 
rate-limiting step is the annihilation of e — h bound pairs, while the binding dynamics is 
much faster. There are different ways of combining UD and BD processes together, and we 
tried some of the simplest possibilities (see below). The model that gave the best results is 
based on the assumption that there is a very rapid initial branching of the charge carrier 
populations into two separate non-communicating channels, respectively decaying with a 
UD and BD law. In other words, we assume that, immediately after the UV excitation, 
two separate charge- carrier populations of electrons and holes are formed, with densities 
Ni(t) and N 2 (t) (hole and electron densities remain equal in each population). We further 
assume that this initial separation has a fixed branching ratio, so that, for any excitation 
energy density U, at time t = we have iVi(O) = a.\U and iV 2 (0) = a 2 U, with a± and a 2 
constants (this is valid only for excitation energies below the saturation level). Next, we 
assume that N\ is governed by a standard UD rate equation, dNi/dt = —Ni/tx, where T\ is 
an intensity-independent characteristic decay time. For N 2 we assume that Eq. (CQ) remains 
valid. Finally we introduce two different fluorescence quantum yields Q\ and Q 2 for the two 
channels. This model leads to the following predicted luminescence decay: 

/ x ~ dNi „ dN 2 tu Co 

'M = -*-3r-*-3r=* e + oW< (2) 

where we have introduced the UD and BD amplitude coefficients 

ci = Q^^/r^Q^U/n 

c 2 = Q 27 iV 2 2 (0) = Q 2l (a 2 U) 2 , (3) 

and the (excitation-energy dependent) characteristic BD time r 2 = l/( r ya 2 U). Expression 
(J2J) (after convolution with the response function) provides an excellent fit to the data, 



as shown for example in Fig. [TJ where the UD and BD terms are also drawn separately. 
Further validation of our model was obtained by requiring that Eq. ([2]) should successfully 
fit the whole data set collected as a function of excitation energy, with a single choice of 
the adjustable parameters (7012) and T\. We note that such parameters determine the BD 
and UD decay times for each given excitation energy U, and are therefore independently 
constrained by each single decay. The excellent fits reported in Fig. [3] demonstrate that 
(7CK2) and n, corresponding in our model to material properties, are indeed independent 
of excitation energy, thus providing a strong confirmation of our model validity. This same 
global fit procedure allowed to rule out alternative models. 

Besides the successful model, we tried the following alternative ones: (i) a pure BD 
equation such as Eq. (jTJ), but acting on unbalanced electron-hole populations (e.g., N e (t) = 
Nh(t) + C), as could arise due to trapping or impurity doping; (ii) a balanced undivided 
population N(t) of electrons and holes obeying the single rate equation dN/dt = —'yN 2 — 
N/t\. The first model cannot fit the single decay. The second can fit the single decay, but 
not the whole set. 

Moreover, we tested a model in which a trimolecular-decay process of the kind dN/ dt = 
—CaN 3 , such as the Auger effect, takes place in the initial stage of the decay, but found that 
it could not fit the observed excitation-energy dependence. Since the Auger effect is entirely 
non-radiative, its presence should also result in a saturated quantum yield versus excitation 
intensity, which in our case starts to appear only above about 1 mJ of excitation energy. 
We may conclude that the role of Auger effect is negligible in our experiments, at least for 
pulse energies below 1 mJ. This result allows us to set the following approximate upper 
limit for the Auger recombination coefficient: C^^IO -33 cm 6 s _1 . This value is two-three 
orders of magnitude smaller than the typical values of the Auger coefficient found in indirect 
semiconductors, such as germanium or silicon. 22|, |23j The reason is probably the larger band 
gap of STO, as the Auger recombination coefficient is known to decrease exponentially with 
increasing gap width. 24] 

In our BD+UD model, the best-fit values of 70:2 were found to range from 5 to 10 
/is _1 cm 3 /J, depending on the sample (corresponding to a BD time r 2 = 300 — 600 ps for 
a UV pulse energy of 1 mJ). Since the ratio a 2 /«i is unknown, we cannot determine the 
value of 7 from this result. However, a plausible order-of- magnitude estimate is obtained by 
assuming a balanced branching ratio otijoLx ~ 1, which yields 7 ~ 10~ n cm 3 s _1 (by a similar 
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FIG. 4: (Color online) Amplitude coefficients c\ and C2, versus excitation pulse energy (proportional 
to the energy density U) of the unimolecular (a) and bimolecular (b) terms in the luminescence 
decay. Square dots are the values extracted from data. Solid lines correspond respectively to a 
linear (a) and quadratic (b) best-fit. Vertical scales of (a) and (b) are the same, but c\ data are 
multiplied by a factor of 30. 

reasoning, we can set the following, well denned, lower limit: 7 > 0.5 x 10~ n cm 3 s _1 ). This 
value is two-three orders of magnitude larger than the bimolecular coefficient for indirect-gap 
semiconductors such as silicon, usually ascribed to phonon-assisted radiative and defect- 



assisted Auger recombination processes. [231] but it is about one- two orders of magnitude 



smaller than the value found in direct band-gap semiconductors. [221 . |25| ] The best-fit values 
of 7~i ranged from 11 to 24 ns. The coefficients C\ and c<i giving the initial amplitude of the 
UD and BD terms in the luminescence were adjusted separately for each excitation energy. 
The resulting best-fit values are shown in Fig. [H where the energy dependence predicted by 
Eqs. ([3]) is seen to be well verified in the non-saturated regime. Besides further confirming 
the validity of our model, this shows that both decay channels are non-saturated in this 
energy range. For higher energies, we observe the onset of saturation in both channels (data 
not shown). The measured yield ratio of the BD to UD channels is of 0.6 — 0.7, with 20% 
sample-to-sample variations but independent of excitation intensity, within our experimental 
uncertainties. 

Finally, we tried to resolve possible spectral differences between our UD and BD processes 
by inserting additional filters in the detection line and checking for variations in the relative 
yields or decay times of UD and BD components. In particular we used a bandpass filter 
for the range 390-480 nm, which therefore selects only the blue section of the spectrum, and 
a long-pass filter with a cutoff wavelength at 495 nm, selecting only the green tail of the 
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spectrum. In both cases, no significant changes were found, indicating that the emission 
spectra of UD and BD processes are largely overlapping. 



IV. DISCUSSION 

Let us now discuss the possible underlying nature of the UD and BD contributions seen 
in our data. We start from the consideration that the small photoluminescence absolute 
yield implies that most electron-hole recombination processes occurring in STO are non- 
radiative, i.e., the radiative transitions that give rise to the luminescence represent only a 
small fraction of the total. This in turn implies that for a given homogeneous population 
of excited electrons and holes in STO, the photo-luminescence signal kinetics is mainly 
controlled by non-radiative decay channels, regardless of the fact that it is actually detected 
through the radiative ones. On the other hand, the spectrum of the photoluminescence 
characterizes only the radiative processes. Therefore, one must be careful when trying to fit 
kinetic and spectral pieces of evidence in a single scenario, as they are likely associated with 
distinct physical processes. 

As we mentioned in the introduction, charge carriers in STO are believed to have a 
polaronic nature, owing to the strong polarizability of the STO ionic lattice. Electron 
carriers in STO are very mobile, especially at low temperatures, as expected from large 
polarons, although there is some evidence that localized small polarons also 



exist in the system. [?], Q, |28| The intrinsic mobility of hole polarons is unknown, as at low 
densities STO holes are usually trapped in mid-gap impurities. Therefore the observed BD 
kinetics is to be probably ascribed to the direct non-radiative recombination of conduction- 
band electrons and valence-band holes, with at least one of the two species remaining mobile 
and delocalized. This recombination process may be either phonon- or defect-assisted. We 



emphasize that the intrinsic BD behavior o 



different from that reported in the past. 



13, 



ie STO luminescence reported here is quite 



18| The latter occurs on a much longer (/xs-to- 



ms) time-scale and at much smaller charge- carrier densities and is usually ascribed to carrier 



29] 



trapping-untrapping processes, dominated by extrinsic effects. 

Focusing now on the radiative processes involved in the BD signal, we must ask ourselves 
why the luminescence spectrum peak is red-shifted by about 0.35 eV with respect to the 
(indirect) band-gap (Stokes' shift) and why the luminescence band extends down to more 
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than 1 eV below the excitation energy. Both energies are far too large to be associated 
with the emission of a single phonon in a phonon-assisted inter-band recombination process 
(the highest-energy phonon in STO has an energy of about 0.1 eV). Midgap impurities 
also cannot be invoked to explain this energy shift, as the number of estimated e — h 
pairs well above 10 20 cm -3 exceeds by orders of magnitude the impurity density, so that 
the impurity emission would be saturated, in contrast to our results. Furthermore, the 
expected sample-to-sample yield fluctuation would be much larger. If small polarons are 
indeed formed in STO, the involved lattice relaxation could instead account for the 0.35 

n 

eV Stokes shift of the blue luminescence peak. |12| However, intrinsic lattice defects such as 
oxygen vacancies or dislocations may also play an important role in the radiative processes, 
either by introducing mid-gap electronic states or by locally enhancing the polarizability of 
the lattice, thus allowing the polaron-like self-localization of charge carriers. 

To account for the UD contribution, we must assume that a fraction of photo-generated 
carriers is rapidly converted into bound pairs of electrons and holes, which then decay 
more slowly with the characteristic lifetime T\ of the bound pair. Due to the indirect 
nature of the gap, e — h pair binding is difficult unless localization takes place. Again, 
it is not understood whether this localization is a purely intrinsic phenomenon or it is 
somehow defect-assisted, although the latter hypothesis seems more plausible, as we find 
significant sample-to-sample fluctuations in the UD component yield and lifetime T\. If 
lattice deformation is instrumental in the localization (both in the perfect crystal or close to 
a defect), the localized bound pair may be identified with the so-called "self-trapped vibronic 
exciton" (STVE). 
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301 ] STVE has been often invoked in the literature to explain 



the strong luminescent emission pea king in the green (at 2.4 eV) seen from STO at low 



temperature (below 40 K) 
criticized. 
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13, 



18 1 , although this interpretation has been recently 
19|, |3jJ In our case, since the spectrum of the UD component seems roughly similar 
to the total spectrum, the STVE should be associated with an emission peaking in the blue, 
as was also recently proposed for the case of electron-induced luminescence in STO.|ll| 

V. CONCLUSIONS 

Summarizing, the time-resolved luminescence decay observed from STO under intense 
pulsed UV excitation is well explained by a model that assumes the presence of two separate 
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decay channels, one associated with a direct recombination of unbound mobile electron- 
hole charge carriers and the other with the recombination of bound electron-hole pairs. 
The exact nature of the microscopic recombination processes taking place in this system 
remains uncertain, although the polaronic character of the charge carriers and the presence of 
intrinsic crystal defects are both likely ingredients of a future understanding. Our experiment 
also demonstrates the possibility of probing this photoinduced high charge- density regime 
without damaging the sample, thus opening the way to several other possible investigations 
of this interesting regime, which may eventually help us clarifying the exact nature of the 
charge carriers occurring in these ever-surprising materials. 
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